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Phyllotaxis on surfaces of constant Gaussian curvature 
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A close packed organization with circular symmetry of a large number of small discs on a plane is 
obtained when the centres of the discs are distributed according to the algorithm of phyllotaxis. We 
study here the distributions obtained on surfaces of constant Gaussian curvatures, positive for the 
sphere or negative for the hyperbolic plane. We examine how the properties of homogeneity, isotropy 
and self-similarity typical of the distribution on the plane, and resulting from the presence of circular 
grain boundaries with quasicrystalline sequences, are affected by the curvature of the bearing surface. 
The quasicrystalline sequences of the grain boundaries appear indeed to be structural invariants, 
but the widths of the grains they separate vary differently with the curvature of the surface. The 
self similarity of the whole organization observed on the plane is therefore lost on the hyperbolic 
plane and the sphere. The evolutions of the local order within the grains show no differences except 
on the equatorial belt of the sphere where the isotropy is decreased owing to the symmetry of this 
finite surface around its equator. 
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INTRODUCTION 

The densest organization of small discs on an infinite plane is obtained when their centres are at the nodes of the 
triangular tiling of the hexagonal lattice. All these discs occupy the same area on the plane and have the same local 
environment of six equally-spaced first neighbours which is reproduced all over the plane according to the laws of 
classical crystallography. This crystalline solution remains valid for a limited assembly of small discs within a finite 
compact domain of the plane under the condition that the borders of the domain be aligned along reticular directions 
of the lattice. If not, as for a circular domain, this solution is no longer valid, the circular symmetry of the domain 
being not accounted for by the finite number of rotational symmetries of a crystal. 

It has been shown that the best packing efficiency found in the case of circular symmetry is that obtained when 
the centres of the discs are regularly placed on the spiral drawn by the algorithm of phyllotaxis [l|-|3| . This solution 
approaches circular symmetry, but the area per disc, the number of first neighbours and their distances are no longer 
constant as they are in the case with hexagonal symmetry. Those differences between crystalline and phyllotactic 
organizations were recently analyzed considering not only the positions of the centres of the discs but also the polygonal 
nature of their Voronoi cells [4[ . Such an analysis shows that the phyllotactic best packing efficiency results from an 
interplay between metric and topological distortions leading to a structure with an inflation-deflation symmetry. We 
study here the distributions of points with their Voronoi cells drawn by the algorithm of phyllotaxis onto surfaces 
of constant Gaussian curvatures, positive for the sphere or negative for the hyperbolic plane. We examine how the 
best packing efficiency typical of the distribution on the plane evolves with the curvature of the underlying surface, 
considering both short-range (the shape of a cell, its area and the distances from its centre to the first neighbour 
points) and long-range characteristics (the repetition law governing the organization of the cells). 

The part of this work concerning plane and spherical phyllotaxis was motivated by our interest for the lateral 
organizations of long biological molecules in dense fibres with a circular section imposed by their surface tension [5[ . 
These organizations can indeed be represented by distributions of points on planar or spherical bases according to 
the fact that the molecules can be parallel or twisted within the fibres. Spherical phyllotaxis was also used earlier 
as a simple ideal approached by vegetal forms |6| , in order to obtain the most homogeneous sampling of points on 
a sphere in numerical integration [7[, to develop climate models of the earth |8| or to estimate the Earth coverage 
of satellite constellations [9(]. We do not know of any application concerning phyllotaxis on surfaces with negative 
Gaussian curvature up to now. 

Since the distributions on the sphere and the hyperbolic plane will be compared to that observed on the plane, we 
recall first the main features of the latter, limiting ourselves to what is needed for this comparison, with demonstrations 
and more details given in |4( . 

PHYLLOTAXIS ON THE EUCLIDEAN PLANE 

The Voronoi cells, used in [4j to describe phyllotaxis on the plane are an efficient tool to analyze the geometry and 
the topology of the structure. 

Pattern 

A phyllotactic distribution of points indexed by s is described by an algorithm such that the position of point s is 
given by its polar coordinates: 

p(s) = a^fs and 9(s) = 2ir\s (1) 

which is the equation of a Fermat spiral called generative spiral hereafter. The parameter a defines the metric scale. 
Sites, indexed by positive integers s, are placed on this spiral, so that the azimuth between two successive points 
varies by 2ir\. The best packing efficiency is obtained for A — 1/r where r is the golden ratio (l + v / 5)/2. The Voronoi 
cell of point s is defined as the region of space nearer to it than to any other point of the pattern. The phyllotactic 
pattern of a set of 3000 points with their Voronoi cells is shown in figure [TJ 



1 In this paper, the first point correspondto s = and the next, to successive integers. Taking s £ Z + 1/2 leads to a more regular 
structure in the core of the phyllotaxis and it appears as a tt/t rotation for large s. 




FIG. 1. A set of n = 3000 points organized on the plane according to the algorithm of phyllotaxis with the golden ratio. Each 
point is surrounded by its polygonal Voronoi cell whose number of sides corresponds to the number of first neighbours around 
the point. Blue, red and green cells are respectively pentagons, hexagons and heptagons. The three spirals joining the first 
neighbours of the points are called parastichies. The white dot marks the origin point s — 0. 



In such a pattern, pentagons and heptagons are topological defects distributed among hexagons. These defects 
appear concentrated in narrow circular rings with constant width, separating large annuli of hexagons whose width 
increases from the core towards the periphery. In the narrow rings, pentagons and heptagons are paired as dipoles 
separated by hexagons whose shape is that of a square with two adjacent corners cut off. The defect rings are 
indeed grain boundaries separating grains of hexagonal cells, each dipole acting as a dislocation that introduces a 
new parastichy (the lines joining first neighbours) to maintain the density as constant as possible. The underlying 
arithmetic of this organization is that of the Fibonacci sequence, f u = f u -i + fu-2 from /o = and f\ = 1, as 
summarized in table |U So cells are organized in concentric blocks, grains separated by grain boundaries, surrounding 
a central core containing about ten points and showing an apparent disorder recalling the structure of confined two- 
dimensions foams. In this paper we are interested in large structures, so that we shall consider what happens outside 
the core: A succession of large grains of hexagonal cells that are concentric circular annuli, bounded and separated by 
circular grain boundaries (f u —i, fu-2, fu—i) made of / u _i hcptagonal cells, f u -2 hexagonal cells and f u -x pentagonal 
cells. 

These grain boundaries serve as natural boundaries for our optimal packing problem. Outwards packing begins with 
the first complete grain boundary (13,8, 13) with 13 heptagons, 8 hexagons and 13 pentagons. The core is bounded 
by the 8 pentagons of the (first) incomplete grain boundary (3,5,8). It has 3 heptagons instead of the full 8. The 
additional pentagons (nearly) fulfil the topological requirement that a tiled circular domain should have a topological 



charge of 6 (i.e. 6 additional pentagons, a sphere having a topological charge 12) [ll|,lF2|. This solves the packing 
efficiency problem: one grain boundary constitutes a perfect circular boundary for the domain into which objects are 
to be packed. 



TABLE I. a) Cell types in the successive rings. The number s is that of each point on the generative spiral, the first neighbors 
of point s have numbers s + Ss where Ss are all Fibonacci numbers except for s = O.The large hexagonal rings or grains are 
labeled by the rank u of the Fibonacci number f u corresponding to the medium value of Ss in a ring, the Ss are also the 
numbers of parastichies of each type in the ring. The narrow rings or grain boundaries are marked by ||. b) Here the core in a) 
has been cleared up by starting the numbering from s — 1 (eliminating cell s — by removing the boundary 0/2) and making 
two neighbour flips to remove the 13 parastichy from the innermost 3 cells. All Ss are Fibonacci numbers [la ]. 



cell type 



number of cells 



s from 



to 



neighbour separations Ss 



pentagon 

hexagon 

heptagon 

hexagon 

pentagon 

hexagon 

hexagon 

heptagon 

hexagon 

pentagon 






1 


1,2,3,4,5 or -1,2,3,5,8 


2 


2 


-2,2,3,5,8,13 


3 


5 


(-3,-2,2) or (-4,-3,-2) or (-5,-3,-2),3,5,8,13 


6 


6 


-5, -0,0, 5, 8, lo 


7 


8 


-5,-o,5,o,lo 


9 


9 


-8,-5,-3,5,8,13 


10 


14 


-8,-5,5,8,13,21 


15 


17 


-13,-8,-5,5,8,13,21 


18 


22 


-13,-8,-5,8,13,21 


23 


30 


-13,-8,8,13,21 



i-1 



hexagon 
heptagon 
hexagon 
pentagon 



fi-i 

fi-2 



-fi, —/i-1, —fi-2, fi-2, j "i-1, fi 

-fi, ~ /i-1, —fi-2, fi-2, /i-1, fi, fi+1 

-fi, — /i-1, —fi-2, /i-1, fi, fi+1 

'fi, —fi-1, fi-1, fi, fi + 1 



b) 



pentagon 

hexagon 

hexagon 

hexagon 

pentagon 

hexagon 

heptagon 

hexagon 

heptagon 



1 


1 


1,2,3,5,8 


2 


3 


-1,1,2,3,5,8 or -2,-1,2,3,5,8 


4 


5 


-3,-2,3,5,8,13 


6 


7 


-5,-0,0,5,0, lo 


8 


8 


-5,-o,5,o,lo 


9 


9 


-8,-5,-3,5,8,13 


10 


10 


-8,-5,-3,5,8,13,21 


11 


15 


-8,-5,5,8,13,21 


16 


17 


-13,-8,-5,5,8,13,21 



Defect rings 



In the middle of hexagonal grains it is easy to identify three parastichies which are spirals running through neigh- 
bouring points. The three curves intersect with angle close to 27r/3, but when the spirals approach the grain boundary 
domain only two remain apparent. The third disappears to restart in an orthogonal direction through the truncated 
square hexagons of the grain boundary. The two apparent parastichies are orthogonal and give the directions of the 
dipolcs. The Voronoi cells on the grain boundaries are nearly squares. The two apparent parastichies are orthogonal 
to the square edges, the third runs along one diagonal only and there is a flip of the diagonal across the hexagons of 
the grain boundary. 

The elementary dipoles (in which heptagon and pentagon are nearest neighbours) are oriented along parastichies 
f u refereed by u in table U but there is a flip of orientation from one grain boundary to the next. Nevertheless all 
dipoles make the same angle arccot(/ u // u _i) ~ arccot(r) in absolute value with the radial direction. This results 
simply from the description of grain boundaries given in appendix A as a strip in a square lattice. 

The perimeter of a ring of dipoles is determined by its number of dipoles and hexagonal cells constituting it. 
Table fl] shows that this number belongs to the Fibonacci sequence so that the ratio of the radii of two successive 
rings of dipoles is approximated by the golden ratio r. Strictly, as shown in appendix A, this ratio is \/ fiu+zl fiu+\ 
converging to r for large Fibonacci numbers. It can be seen on figure [1] that the dipoles occur as singletons, or 
in pairs, distributed on the ring according to an approximant of a quasicrystalline sequence |19| . as explained in 
appendix A. When moving from one ring to the next, away from the core, the evolution of the sequence is determined 



FIG. 2. (a) Variation of the area of Voronoi cells for a pattern of n — 3000 cells. This is obtained with the scaling parameter 
a = 1 in equation [1] leading to an area close to it. This area tends toward it, but with rapid variations in grain boundaries. 
Inserted is an histogram of the area of cells, notice that the scale is very inflated, (b) First neighbour distances between points 
s and s + Ss. Blue, for the interval 5s equal to the smaller positive Fibonacci number in the list of table 1. Green, for the 
next interval, red for the third positive interval and purple for the last one (occurring only if the Voronoi cell is an heptagon). 
So green, blue and red correspond to distance along the three visible parastichies. Each continuous curve corresponds to a 
given Fibonacci number that appears in different rings. For instance, /n = 55 appears between s — 101 and 2254 leading to 
a continuous curve which is successively purple, red, green and blue. Lower and upper crossings of two curves corresponds 
to a grain boundary, other crossings are in the middle of a hexagonal grain. An histogram of all distances is inserted. The 
hook for distances close to 1.9 and 2.3 is an artefact due to the finite number of points. The distances lies in the range 
[(2n/V5) 1/2 ,(2^% 




1000 1500 1000 2500 3000 500 1000 1500 2000 2500 



300(1 



by the Fibonacci inflation/deflation rule where a singleton becomes a pair and a pair becomes a pair plus a singleton. 



Metric properties 



The area of a disk of radius p which contains s points is irp 2 = ira 2 s so that the average area per point has the 
value 7ra 2 . The Voronoi decomposition is a tessellation which attributes a specific area to each point. But the areas of 
Voronoi cells are not all identical as shown in figure [2^, even if fluctuations are small. The area of a Voronoi cell falls 
rapidly each time a ring of dipoles is crossed as s increases and it tends towards ira 2 for large s. The fluctuations are 
small immediately outside the core of the pattern. The histogram describes the level of homogeneity attained with 
phyllotaxis. The standard deviations of similar histograms arc 8 — 0.01589 for n = 6000, S = 0.02246 for n — 3000 
and 5 = 0.03171 for n = 1500 (assuming a = 1), rapidly decreasing with size n. 

The behaviour of the distances between first neighbour points is shown on figure (2b- In grain boundaries where 
Voronoi cells are slightly deformed squares, the cell area is approximatively ira 2 so that the distance between points, 
along a square edge, is y/n ~ 1.772 (with a = 1) and along a diagonal, is \/2tt ~ 2.506. In the grains where Voronoi 

cells are more clearly hexagons two distances are close to « / ■% ~ 2.053. The distance between the two points defined 

by s and s + f u as a function of s is obtained using the relation (see [4J, |l4[): 



d u (s) = f u (— -+■ 



1 , s{-2Kf u _ 1 + 2vT- l f u f 



t-i 

J U 



)1/2) 



(2) 



this relational fits the numerical values given on figure [2}d. It can be checked that minimal values are 



V^|i-^l (3) 

all close to and converging toward J^j? — 1.67. All distances between first neighbours are in the range [\/-y?, \/27r], 

whatever will be s and u. This property, specific to the choice of A = 1/t is very important to ensure the best 
uniformity in all orientations. The reason is that t is approximated by successive truncations of its continuous 



fraction (having only 1) which are the ratio of two successive Fibonacci numbers and so converging smoothly [16l.ll7[. 
It is this property which leads to have minimal values for distances given by the equation [3] converging toward a 
finite value. Other ratios (with a tail of their continuous fraction expansion containing other integers) introduce some 
distances between neighbours decreasing rapidly with the Ss separation between neighbours. 

PHYLLOTAXIS ON THE HYPERBOLIC PLANE 
The Poincare disc 

The Poincare disc model is a simple way to represent the hyperbolic plane as presented in appendix B. A point 
in the hyperbolic plane is defined by polar coordinate {p, i?) and is represented on the Poincare disc by (r, ■&) where 
lengthes are obtained with a given metric da 2 — 4(dr 2 + r 2 dd 2 )/(l ~ r 2 ) 2 . So r = 1 represents points at infinity. 
This metric has the form of an Euclidean metric in polar coordinates divided by a function of r only: it is locally an 
Euclidean metric. 

Pattern 

As for phyllotaxis in the Euclidean plane, a point is defined by an integer s on a spiral and the pattern obtained 
is shown in figure [3] The number of points in an hyperbolic cap of radius p(s) is proportional to s. We choose the 
p(s) function so that the number of points enclosed in the domain is ira 2 s, as it is for the Euclidean plane phyllotaxis. 
Using the equation [9] given for the area enclosed by a circle in appendix B, we write 2ira 2 s/2 = 2-7r(cosh<p — 1)> with 
R = 1 to have the Gaussian curvature equal to n = —1, then ip = p{s). The true radius in the hyperbolic plane is 

p(s) = cosh -1 I %p + 1 1 . This radius is a length measured on the hyperbolic plane, but on the Poincare disc, The 

Euclidean distance from the origin to the representation of a point r = tanh(<^/2) is given by 

r(s) = tanh ( - cosh -1 ( — + 1 ) ) . (4) 

The factor a 2 / 2 relating s and p(s) allows to have the radius in the hyperbolic plane p(s) ~ ciy/s for small s 
as in flat phyllotaxis. The spiral equation in the hyperbolic plane is (p(s) cos(2tt Xs) , p(s) sm(2n Xs)) where A is a 
parameter exactly similar to the A parameter introduced in plane phyllotaxis. Similarly the true hyperbolic phyllotaxis 
corresponds to the choice A = 1/t. As shown in 3.4, the p(s) behaviour and the choice of A = 1/t ensure the best 
homogeneity and isotropy for the distribution of points. 

Rings of dipoles 

A Voronoi decomposition of the phyllotaxis mapped on the Poincare disc can be obtained as if the structure in this 
disc was Euclidean. This is justified by the fact that the metric is locally Euclidean. If two points are neighbours with 
the hyperbolic metric they are also neighbours using the local Euclidean metric on the Poincare map lj. On figure El 



This relation suppose that Ss is infinitesimal compared to s. An estimation of the accuracy can be obtained comparing the distances 
computed between points s and s + Ss with that between s + Ss and s. They must be equal but computed values are slightly different. 
An average between both is very accurate. 

More precisely, this is also related to the fact that two triangular Delaunay decompositions obtained using the Euclidean metric of 
the plane representation or the hyperbolic metric are formed with the same set of triangles. These sets of triangles are defined by the 
circumcirclcs of triangles selected so that they do not enclosed any points of the phyllotaxis. Because circles in the hyperbolic plane 
are mapped onto circles in the Poincare representation the two selected sets are identical. Consequently neighbouring relations are the 
same with the two metrics. 



FIG. 3. Hyperbolic phyllotaxis represented inside the Poincare disc. The pattern is obtained using equation[4]with 3000 points 
and a = 1/20. Blue, red and green cells are respectively pentagons, hexagons and heptagons. All area and distances are scaled 
so that the mean area per points is tt. With this choice the curvature k = —1/R 2 is given by R = 1/a, with the parameter 
a — 1/20. By effect of the Poincare representation size of cells seems to decrease going toward the limit circle. 




FIG. 4. Area and distances on the hyperbolic plane. All area and distances are scaled so that the mean area per points is tt. 
With this choice the curvature n = —1/R 2 is given by R — 1/a, with the parameter a — 1/40. (a) Area of Voronoi cells for the 
hyperbolic phyllotaxis. The area are close to the mean value tt, strong variations appear near rings of defects. Nevertheless 
these strong variations decrease as the area tends toward n for large s. (b) Distance between point s and its first neighbours 
s + Ss with positive Ss. Green, blue and red correspond to distance along the three visible parastichies. Distances are confined 

between the limits y/2n and J 2-k j\[Z. 
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Voronoi cell edges are straight segments; on an exact Voronoi decomposition they must be arc of circles representing 
geodesic lines of the Poincare disc but as these segments are very short, the approximation is clearly good. 

As in plane phyllotaxis it is possible to define hexagonal grains separated by grain boundaries formed with pentagon- 
heptagon dipoles. It is the change in the list of separations Ss = f u -i, fu, fu+i between neighbours in a grain and 
in the following one which govern a grain boundary. So the description given in appendix A remains the same, still 
related to Fibonacci numbers. Grain boundaries are indeed the same in the three examples of phyllotaxis. 



Metric properties in hyperbolic phyllotaxis 

The p(s) function has been chosen so that the area enclosed in a circle of radius p is ira s on a hyperbolic plane of 
K = — 1 Gaussian curvature. On figure H^, the evolution of the area per cell is shown scaled in order to have the mean 
area equal to 7r as in the plane case; so the curvature is k = —1/R 2 with R = 1/a. Area fluctuations are stronger, 
with a smaller damping than in the plane case. 

It is possible to get distances between first neighbour points from the coordinates; this is presented on figure 0p 
but it is also possible to have an accurate analytical relation allowing to follow behaviour of distances with s. On a 
plane representation, the distance, measured on the plane, separating two close points whose indices are s and s + Ss 
with Ss = f u is given with a good accuracy by: 

{d u {s)f = f u [r'{s) + { lu /f u fr{sf], (5) 

with 7„ = (— 27r/„_i + 2iTT~ 1 f u ) and r(s) given by equation |4] Equation [2] given for a plane phyllotaxis is derived 
from this equation (see the foot note in 2.3 on the accuracy of this relation), but to have the true hyperbolic distances 
requires to take account of a correcting factor to the metric. This factor is 1/(1 — r(s) 2 ). Then the distances are given 
by: 

, , v , _1 , -a 2 , r( S ) 2 (-2^/. lt _ 1 +27TT- 1 / 1[ ) 2 u/2 

duis) = fu r^w { Mi + sa 2 r + Ji } ■ (6) 

This relation fits the curved given on figure 0p after normalization by a scaling factor 1/a in order to have mean area 
of cells equal it. 

As the sequence of rings of dipoles are the same whatever the curvature as developed in appendix A, the effect of the 
negative curvature is to decrease the width of hexagonal grains which tends to be a constant for large s. Effectively, 
the radius of a grain boundary is sinh^u = y\^2u+i/7r/2 which leads for large ip u to ip u ~ ln(\/f2u+i/^)- Two 
successive grain boundaries defined by u — 1 and u give a width ip u — ip u -i ~ ln(\//2u+i//2u-i) very close to the 
constant In r. 

Like in the case of the plane, the maximum for first neighbour distances is still y2n ~ 2.506 corresponding, to 
cells close to squares with a distance corresponding to the diagonal. The minimal values of d u (s) (equation [6|) are 

V^fu\/\- — f ± \ very close to and converging toward . /-^= ~ 1.67. Distances are confined in the same domain as 

for the plane phyllotaxis, a consequence of the choice of the parameter A = 1/r. 

PHYLLOTAXIS ON THE SPHERE 
Pattern 

The algorithm used to build the phyllotactic configurations shown in figure [5] is such that the position of point s is 
given by the spherical coordinates (p,cj>,0). These coordinates are p = R the sphere radius, 8 = 2irAs, the azimuthal 
angle and <j> = arcsin(s'/i>) — n/2, the polar angle. The total number of points on the sphere is n = 2v + 1. With 
this choice, the integer s' goes from —v to v with s' = —v or s' = v for points on poles, and s' = for a point on the 
equator. For comparison with plane phyllotaxis it could be helpful to have s — on the north pole using s = s' + v. 

Spherical phyllotaxis is obtained by mapping of points on a finite cylinder which is tangent to the sphere along 
the equator, so which has the same radius R and whose finite hight is 2R, the area of this finite cylinder is 4irR 2 as 
that of the sphere. The points on the cylinder are located on a perfect helix defined by 9 — 2ir\s and a pitch related 
to the number of points on the sphere. They can be consider as a perfect crystal without defect wrapped on the 
cylinder, so with a constant density of points. This is an example of phyllotaxis on a cylinder as described in [la. Hal- 
With A = 1/r, a point s has six neighbours at position s ± Ss with Ss equal to three successive Fibonacci numbers, 
depending on the pitch of the helix drawn on the cylinder. These points are mapped on the sphere orthogonally to 
the polar axis and so the area per points on the sphere is the same as on the cylinder. This results from the property 
of the area of a spherical zone enclosed on a sphere of radius R, between two parallel planes at distance h to be the 
same as the area of a finite part of a cylinder of radius R enclosed between the two same parallel planes orthogonal 
to the cylinder axis. This justify the choice of the polar angle on the sphere <j> = arcsin(s'/^) — 7r/2 deduced from the 
cylindrical coordinate (R,9, Rsm(cf) + tt/2)). But the projection on the sphere, even if it conserves area, introduces 
inhomogeneous shearing in the structure. This shearing, at constant area, is given by a compression factor cos(0) 



FIG. 5. Phyllotaxis on a sphere with 4001 points. Blue, red and green cells are respectively pentagons, hexagons and heptagons; 
the white point is a pole. Ring of dipoles appear clearly. Near the equator there is a large hexagonal domain even if cells look 
like squares, in fact they have six edges two of which being very small. The structure has a chiral symmetry around the axis 
joining the two poles. 




changing the equator length into a parallel length (at <f>), and an expansion factor l/cos(0) along meridians. This 
shearing could change the incidence relations between neighbours so first neighbours are not necessarily the same 
as on the "crystalline" cylinder and defects appear. The area associated to cells is close to be constant as in plane 
phyllotaxis. The choice of R is arbitrary but in order to have cell area close to it as in plane phyllotaxis R = \J2v + 1/2. 
A Voronoi decomposition of the set of points on the sphere reveals the phyllotaxis and how defects appear Q As 
in plane phyllotaxis the Voronoi cells are hexagons, pentagons or heptagons, with most of the Voronoi cells being 
hexagons. We call the large domains with only hexagonal Voronoi cells in hexagonal grains, even if they are not 
regular, and the circular distribution of heptagonal, hexagonal and pentagonal cells are grain boundaries. 

Rings of dipoles 

The number of rings of dipoles depends on the number n of points on the sphere, related to its radius R = y/n/2. 
Figure [5] shows how new dipoles appear increasing n. The evolution of rings of dipoles is related to the most 
important fact that the distributions of dipoles and hexagonal cells along these rings follow the one dimensional 
sequences approximant of quasicrystals as developed in appendix A. The number of cells in these rings is 2f u + /«-i, 
where f u is the number of dipoles, and f u -i, the number of hexagons. The length of circles defining rings of defects 



4 The Voronoi decomposition of a set of points on a sphere has been done using a 3D Voronoi analysis considering the set of points on the 
sphere to which is added the center of the sphere. Then Dclaunay tetrahedra are obtained. They are formed with Dclaunay triangles 
on the sphere surface completed by the center of the sphere. The Voronoi decomposition on the surface is deduced from this Dclaunay 
triangulation of the surface. 
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are the length of the strip in a square lattice, so this length is the same in plane and spherical phyllotaxis. It is 

\/(fu + fu+i) n or \//2u+i7!" if we suppose an area of cells given by n and then a cell edge of ^/tt when cells are 

squares. Here and in the following we set the sphere radius R = \]2v + 1/2 in order to have average cell area equal 
to 7r as in plane phyllotaxis. This is then used to estimate the number of cells in hexagonal grains. 

On a sphere of radius R = \/2v + 1/2 containing n = 2v + 1 sites, a circle defined by the polar angle <j) border a 
spherical cap containing v(\ — coscj)) sites. If the circle corresponds to a grain boundary of length \J f2u+i^ this leads 

to a polar angle <fi given by sin0 = \/ (iuVx\ ' ^ e estimate the enclose number of points inside the spherical cap 
bordered by the grain boundary to be: v(\ — i 1 — ,„ Ji °t! ). This is an estimation, in order to get size of hexagonal 



{2u+1)tv> 

grains we have to take account of the width of the grain boundary and suppose that the estimation corresponds to 
the medium of the f u —i hexagonal cells in the grain boundary. An estimation for grain boundaries bounds, sometime 
shifted from 1 is: 



[(5-/ u _ l )/2 + Kl-Jl- ,/ 2 " + 1 \ )J,L(/„-i +3)/2 + i/(l- \h~ J 2 T\ )J - (7) 

y (2i/ + ij7r y (2i/ + ij7r 

Increasing the number n of points on the sphere, and so decreasing the curvature, new rings of dipoles (grain 
boundaries) appear on the equator <f> = tt/2. These is presented on figure [5] in the range n 6 [1331,1351]. New 
grain boundaries appear by pair: one in each hemisphere. As they have a finite width just after the threshold of 
appearance the two grain boundaries arc interpenetrating; this leads to some cells which could have four edges. Setting 
sin (j) = 1 gives the the threshold for the number of points on the sphere for which a new grain boundary appears on 
the equator. The first set of these values are n = (1, 2, 4, 11, 28, 74, 194, 508, 1331, 3484, 9122, 23881, . . .). When a new 
grain boundary, related to the Fibonacci number /2U+1, appears the previous one, related to the Fibonacci number 

jW-ii enters a polar cap of angle <ft given by sin0 = <J ^ ~ 1 or sine/) = \ f~ l ■ It results that cf> ~ arcsinr -1 about 

0.666 rd. This value of <j> does not depend on n as soon the approximation \/ f 1 — T ~ X is good. 

The same evolution repeats itself quite regularly at each appearance of a new grain boundary as shown on figure [7] 
Using a semi-logarithmic horizontal scale, the curves can be translated from one to the next. 

Metric properties on the sphere 

Figure [8] shows that, each time the number n of point on the sphere is translated from one value to the next 
along the logarithmic scale of figure the distances between first neighbour points are reproduced while a new grain 
boundary is introduced. The local order varies within the same limits whatever the number of points. Far away from 
grain boundaries the area are close to the mean value n corresponding the the choice R = \fnj1. There are strong 
variations crossing grain boundaries near the poles which decrease rapidly. 

In hyperbolic case we have evaluated distances in the Poincare representation using a given metric. It is convenient 
to use a similar method for the spherical space using a representation of the sphere as described in appendix C, 
which is no more than a stereographic projection of the sphere on a tangent plane. On this plane representation of a 
spherical phyllotaxis the distance separating two close points whose indices are s and s + Ss with Ss = f u is always 
given by the equation [SJ but where r(s) is given by the equation [TUJ With equation [5] this distance is evaluated with 
the plane metric, so to have the true distances on the sphere needs to take account of a correcting factor to the metric 
which is 1/(1 + r(s) 2 ) as given in appendix C by equation [TTJ Then the distances are given by: 

-.(.) = /.'( .. -°' . 2]3 + ""'<- */-- + »"-/.)> (8) 

1 + rysf 4s(— 1 + sar)- 3 /„ 

This analytical distances fit the distances presented on the figure [5] obtained from coordinates of the spherical phyl- 
lotaxis taking account of a scaling factor 1/a in order to have mean area of cells equal tt. The important conclusion 
is that these distances are again confined in the same domain found for the plane phyllotaxis. Like in the case of 
the plane, the maximum for distance is still v2tt ~ 2.506 corresponding, to cells close to squares with a distance 

corresponding to the diagonal. It can be checked that minimal values of d u (s) (equation [5]) are V^fu\/\- — f- L \ 

very close to and converging toward \j -js — 1-67 like for the plane phyllotaxis. As in the plane and hyperbolic cases, 
it is this confinement of distances which is the signature of the best uniformity in all orientations. 
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FIG. 6. Evolution of rings of defects near the equator increasing the number n — 2v + 1 of points on the sphere. Colours 
of Voronoi cells are red for hexagons, green for heptagons and blue for pentagons. Appearing on the equator they are two 
interpenetrating rings. Then some cells having only four sides appear in yellow. 
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FIG. 7. a) Evolution of the polar angle cj> of grain boundaries with the number n of points on spheres (logarithmic horizontal 
scale). Each curve corresponds to grain boundaries defined by u — 5,6,7,8,9,10,11,12. For a given number n (or a given 
radius R — y/n/2) there is a grain boundary near the equator if there is a curve such that <f> ~ n/2. The two black points refer 
to the two examples of spheres whose cells area and distances are given on figure [S] b) Evolution, with the Gaussian curvature, 
of the arc length which is the radius of grain boundaries (for instance Rxj> in spherical case). There is a limit to these curves in 
spherical case corresponding to the appearance of grain boundaries at the equator. 
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CONCLUSION 



The spiral organization generated by the algorithm of phyllotaxis with the golden ratio ensures the best packing 
efficiency of points on a plane in a situation of circular symmetry |3(. We examined here how this simple schema 
could be affected when the bearing surface presents a positive or negative Gaussian curvature, limiting ourselves to 
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FIG. 8. (Left) Area of Voronoi cells for the spherical phyllotaxis. All area and distances are scaled so that the mean area 
per points is n. Two examples are given with n = 9301 and n = 1351 points on the sphere. (Right) For the same spheres, 
distance between points s and its first neighbours s + Ss with positive Ss in the north hemisphere and negative 8s in the south. 

Green, blue and red correspond to distance along these three visible parastichies. Distances are in the limits \^2ty and \/2n/-\/E. 
Curves for the two examples n = 9301 and 1351 look very similar, mainly for distances: this is because the n values are chosen 
to have grain boundaries at the same angular positions as given on figure [7] 
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curvature radii much larger than the mean distance between points. 

Drawings of figures Q] [3] and [5] show that the phyllotactic patterns always present the same appearance. They are 
all made of a core, without any obvious discernable order, surrounded by an ordered series of alternate concentric 
rings. The cores are the same, they occupy a limited region around the centres of the patterns where the curved 
surfaces can be assimilated to their tangent plane. The alternate concentric rings are respectively large rings, or 
grains, containing points with six first neighbors and narrow rings, or grain boundaries, containing also points with 
five and seven first neighbors associated in dipolcs, or dislocations. Those topological defects are needed to maintain 
the density as constant as possible in this situation of circular symmetry. 

A detailed analysis of grain boundaries makes apparent the fact that they can be found identical to themselves 
on the three kinds of surfaces. The number, distribution and orientation of dipoles on each of them, as well as 
their evolution from one to the next, do not depend on the curvature of the surface, their perimeters just follow 
the Fibonacci series as shown in table I. This directly proceeds from the mathematical structure of the algorithm of 
phyllotaxis with the golden ratio as developed in appendix A. Figure [9] shows indeed that the dipoles are organized 
along a grain boundary according to a quasicrystallinc sequence whose evolution from one grain boundary to the 
next is driven by a well defined inflation/deflation rule. The distributions of the whole set of points on surfaces of 
different curvatures must respect this most remarkable structural invariance of the grain boundaries as developed in 
the following paragraphs. 

As the perimeters of the grain boundaries are constrained to follow the Fibonacci series, the distance between them 
on the surfaces, or the widths of the grains they enclose, vary differently according to the sign of the curvature. If we 
call P the perimeter and R the radius of curvature, the radii of the grain boundaries indeed varies as P/2ir on the 
plane, i?arcsin(P/27ri?) on the sphere and i?argsh(P/27ri?) on the hyperbolic plane and the width of a given grain 
decreases as the curvature goes from a positive to a negative value. This results in the more or less rapid dampening 
of the area per cell around ir visible on figures [2] and 0] 

As the dipoles on the grain boundaries have well-defined alternate orientations relative to the radius vector with an 
angle close to 7r/2 — arctan(r), the anchoring conditions of the parastichies at the two limits of any grain are the same 
whatever the curvature of the surface and the width of the grain. This constraint manifests itself by the fact that the 
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distances between first neighbors oscillate, along more or less extended horizontal scales, between the same limits and 
cross each other at the same level on the whole plane and hyperbolic plane as shown by figures [2] and [4] However, 
on the sphere, such an oscillating behaviour is observed on a limited polar cap only and not around the equator, as 
shown on figure [3] This holds to the fact that, when two new grain boundaries merge along the equator, their dipoles 
are parallel and stay parallel as they move towards their pole as the sphere grows. The anchoring conditions of the 
parastichies in the grain surrounding the equator do not alternate as they do in other grains and the shape of its 
Voronoi cells evolves differently than in normal grains. Figures [5] and [6] show that cells in the equatorial belt may 
keep a shape close to that of a square rather than becoming hexagonal as in normal grains. 

Finally, owing to the inflation/deflation rule determining the quasi crystalline sequences of the grain boundaries 
they form by themselves a self similar set which is scale invariant as the characteristic distance is changed by a factor 
which is an approximant of r™. When on the plane, this self similarity is transferred onto the pattern as, assimilating 
grain boundaries with ideal circles, their perimeter vary as successive approximants of r and the area of the grains 
vary as r 2 ™ . This is however not valid on the sphere or the hyperbolic plane where the area of the grains do not follow 
a variation in r 2 ™. 



APPENDIX: GRAIN BOUNDARIES AND INFLATION-DEFLATION SYMMETRY 



The way dislocation dipoles are organized along circles is strongly related to ID quasicrystals [19(. A Fibonacci 
ID quasicrystal can be obtained using an inflation-deflation rule iteratively applied to a sequence of long and short 
segments. This inflation-deflation rule is L — >• L + S and S —>■ L. Starting simply from a short segment this specific 
rule gives a quasicrystal after an infinite number of iterations, or with a given number of iterations, a finite structure 
with a number of short and long segments given by two successive Fibonacci numbers. Consider now dipoles formed 
by an heptagon and a pentagon in contact along grain boundaries. There are isolated dipoles, singletons and pair of 
close dipoles. The application of the rule changes a singleton into a pair of dipoles and changes pair of dipoles into a 
pair and a singleton, and transforms one ring of dipoles into the next one. So there is an inflation-deflation symmetry 
associated with radial change relating defects in this structure. This can be checked counting the number of pair of 
dipoles or of isolated dipoles which are successive Fibonacci numbers, on circles of defects. 

One dimension quasicrystals can be obtained by the cut and projection method, selecting the points of a square 
lattice falling inside a strip defined by translation of a cell of the lattice along a straight line of slope r. Approximants 
are obtained when the slope is a convergent of r given by the ratio of two successive Fibonacci numbers. The figure 
[9] which describes how a grain boundary can be derived from a selection of points in a square lattices is similar to the 
cut and projection method applied to approximants. 

The description of grain boundaries as a strip of squares, allows to estimate their lengths. We suppose, as it is 
done in the text, a square edge length d = \JH for square of ir area. The length of the grain with / u+ i dipoles is the 
module of the vector (/„, f u +i) which is L = v/^(/u + fu+i) 1 ^ 2 or simply (f2u+i^) 1 ^ 2 ■ Folding the strip into a ring, 
this length is the perimeter of a circle, then it is possible to estimate the number of points of the phyllotactic pattern 
in the domain enclosed by this circle. 

Dipoles which are oriented along parastichies (refereed by u in table HJ are close to make the same angle with the 
radial direction (in absolute value as their orientations alternate) . Considering the description of grain boundaries by 
a strip in a square lattice (figure ^ , dipoles are represented along square edges and make a constant angle with the 
large strip side related to the slope of the strip f u /fu-i- The angle with the normal to the strip is arccot(/ M // u _i) ~ 
arccot(r) ~ 0.5535 rd. This is also the angle of a dipole with the radial direction for a strip refolded into the ring of 
a grain boundary if we consider that resulting distortions are small. In fact this is exactly the angle with the radial 
direction done by the parastichies u through the medium point of the grain boundary, whatever is the curvature. 

There are two correlated properties of inflation-deflation symmetry related to grain boundaries. The first one is the 
organization of dipole on a grain boundary which is like an approximant. The other aspect, developed in [J] relates 
a grain boundary to the next one. There is a symmetry mixing spiral symmetry and inflation-deflation symmetry 
associating the two grain. Nevertheless such symmetry needs the scale invariance of the plane and cannot be considered 
in curved geometries. In plane geometry the ratio of the radius of two successive grain boundaries is \/ Hu+zl fiu+\ 
converging toward r for large Fibonacci numbers. This is a consequence of the inflation-deflation symmetry in plane 
geometry The structure of a grain boundary is only related to its rank u through Fibonacci number f u and does not 
depend on the space curvature. It is the same for spherical, plane and hyperbolic phyllotaxis. 
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APPENDIX: THE POINCARE DISC REPRESENTATION OF HYPERBOLIC PLANE 

The Poincare disc model, also known as the conformal disc model is a simple way to represent the hyperbolic 
plane. It is a mapping of the whole hyperbolic plane on the interior of a circle, so this limit circle represent infinity. 
It is a conformal mapping respec ting angles between geodesic lines which are represented by arcs of circles that are 
orthogonal to the limit circle [20|, |21| . 

A method to introduce the hyperbolic plane and its representation, is to refer to spherical geometry changing the 
Gaussian curvature k = 1/R 2 into — 1/R 2 , where R is the radius of curvature of the surface. By analogy with the 
equation of a sphere embedded in the Euclidean 3D space, it would be tempting to write x 2 + y 2 + z 2 = —R 2 with an 
imaginary radius, but the hyperbolic plane can not be entirely embedded in K 3 . So we need a representation of the 
hyperbolic plane by a surface with a given metric. A good example is given by the Minkowski plane. Consider a 3D 
vector y = (2/0,2/1,2/2) in a 3D space with a metric such that the squared modulus of the vector is | y | 2 = — 2/o + 2/i +2/i- 
The equation — y 2 + 2/1+2/2 = —R 2 defines a two sheets hyperboloid surface. A single sheet, with the appropriate 
metric, is a model of a surface, with Minkowski metric, whose Gaussian curvature is —1/R 2 . Geodesies of this surface 
are intersection lines between the hyperboloid and planes through the origin. 

Using polar coordinates to describe the y vector we have 2/0 = R cosh ip, y\ = i£ sinh </> cos $, 2/2 = R sinh ip sin •& for 
points on a sheet. The Minkowski metric da 2 = —dy 2 +dy 2 +dy 2 is with polar coordinates da 2 = i? 2 (d</? 2 +sinh tpd-& 2 ). 
By comparison with spherical geometry it is interesting to note that if we change variables into complex variables 
(yo = rjo, yi = irji, 2/2 = ir]2) and (p = i(j> the metric look like the usual metric of a sphere da 2 = R 2 {d(f) 2 + sin (j>d9 2 ). 
In spherical geometry, the parameter R is the radius of the sphere. In hyperbolic geometry R can not be consider as 
a radius of curvature, it is only the Gaussian curvature k = —1/R 2 which is an intrinsical property of the surface. 

The Poincare disc representation is obtained by stereographic projection of one hyperboloid sheet unto the plane 
(2/1,2/2) with a projection pole (2/0 = — 1,2/1 = 0, 2/2 = 0). We consider polar coordinates (r, 1?) in the plane (2/1,2/2) 
(that is in the Poincare disc). The stereographic projection give r = Rsmh(p/(1 + 2/0) or r = Rta,nh(ip/2). Using 
hyperbolic trigonometric identities the metric is da 2 = 4R(dr 2 + r 2 d-d 2 )/(l — r 2 ) 2 . This metric has the form of an 
Euclidean metric in polar coordinates divide by a function of r only: it is locally an Euclidean metric, a proof that 
the representation is conform. 

To obtain the equation of the spiral allowing to build a phyllotaxis on the hyperbolic plane, the calculation of the 
area enclosed in a circle of a given radius R(p p is needed. This area is obtained using this metric. The perimeter of 
the circle centered at the origin is 2irRsmhip and so the area enclosed by this circle is JT P 2irR 2 sinh ipdtp which is: 

27ri? 2 (coshv3 p -l). (9) 

APPENDIX: A PLANE REPRESENTATION OF THE SPHERICAL PHYLLOTAXIS 

There is an other way to describe the spherical phyllotaxis by analogy to the hyperbolic one, simply changing 
hyperbolic trigonometry into spherical trigonometry. This give a representation on the plane of a spherical phyllotaxis, 
with an adjusted metric as shown in figure [TU1 In fact it is a stereographic projection of the tiled sphere onto a tangent 
plane through north pole from a projection point on the south pole. 

The radial position for a point is given on the plane by 

r(s) =tan( -cos" 1 (l -a 2 s/2) j . (10) 

This radial value is the radial coordinate on the plane and the true spherical distance from the origin to the point 
depends on the spherical metric: 

da 2 =A(dr 2 +r 2 d9 2 )/(l + r 2 ) 2 . (11) 

This is a metric which has the form of an Euclidean metric in polar coordinates divided by a function of r only. 
This is practically useful to calculate distances: they are calculated as Euclidean distances on the representation, 
then corrected using the r function. This is a good approximation if distances remain small compare to the radius of 
curvature. This locally Euclidean metric is related to the fact that the stereographic projection on the plane of the 
spherical tiling is conform. The azimuthal position is the same as in the plane case: 2irs\ with A = 1/r. 

Comparing the Euclidean plane case to the plane representation of curved spaces like the Poincare disc in the 
hyperbolic case or the stereographic map of the sphere, the main difference is the behaviour of the radial part 
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r(s) of the spiral equation. Evidcntally if we consider the usual metric of the plane, and not the given metric of 
representations, the size of cells increases in the stereographic map and decreases on the Poincare disc. It seems 
possible to generalize this observation to phyllotaxis where the size of cells is not constant. Any plane phyllotaxis 
build using a spiral equation in polar coordinates (r(s),9(s)) with 9(s) = 2-ks/t would have always the same grain 
boundaries. For instance this apply to composed flowers like sunflower where the size of florets increase from the core 
to the periphery due to growing, very similar to the plane representation of the spherical phyllotaxis but with the 
plane metric. 
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FIG. 9. Strip cut in a square lattice. This strip is a good approximation of the ring of defects containing 55 points. Color 
of the points correspond to the type of their Voronoi cells: hexagon (red), heptagon (green) and pentagon (blue). The first 
point in green is down, then points are numbered, increasing by 21 going up or decreasing by 34 going right. This strip can be 
divided into three strips of heptagons, hexagons and pentagons. 
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FIG. 10. Stereographic projection of a spherical phyllotaxis (part). This is given using equation (CI) with 2000 represented 
points and a — 1/40 so the total number of points on the sphere is 6400. By effect of the projection size of cells seem to increase 
going outward, diverging for n — 6400. 




